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Abstract 

Using special linear combinations of finite energy sum rules which 
minimize the contribution of the unknown continuum spectral func- 
tion, we compute the decay constants of the pseudoscalar mesons B 
and B s . In the computation, we employ the recent three loop cal- 
culation of the pseudoscalar two-point function expanded in powers 
of the running bottom quark mass. The sum rules show remarkable 
stability over a wide range of the upper limit of the finite energy 
integration. We obtain the following results for the pseudoscalar de- 
cay constants: f B = 178 ± 14 MeV and f Bs = 200 ± 14 MeV. The 
results are somewhat lower than recent predictions based on Borel 
transform, lattice computations or HQET. Our sum rule approach 
of exploiting QCD quark hadron duality differs significantly from 
the usual ones, and we believe that the errors due to theoretical 
uncertainties are smaller. 
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1 Introduction 



The decay constant of a pseudoscalar meson B q consisting of a heavy 6-quark 
and a light quark q, with q = u,d, s, is defined through the matrix element of 
the pseudoscalar current: 

<n\(m b + m q )qi Tj b)(0)\B q > = f Bq M 2 Bq . 

These decay constants are of great phenomenological interest since they enter 
in the input to non-leptonic £?-decays, in the hadronic matrix elements of B — B 
mixing and in the extraction of \V c b\ from the leptonic decay widths of B-Mesons. 
Knowledge of the decay constants allows to estimate the so-called hadronic B 
parameter which is directly related to the deviation of the vacuum saturation 
hypothesis. The decay constants are therefore of central interest to the ongoing 
experiments carried out in B- factories. Unfortunately these matrix elements 
could not be measured directly so far, so that we have to rely on theoretical 
calculations. As the calculations must be non-perturbative there are essentially 
two approaches, QCD sum rules and lattice simulations. 

The method of the sum rules has been successfully applied since the pioneer- 
ing work of Shifman, Vainshtcin and Zacharov 1 to calculate various low energy 
parameters in QCD. Particular sum rules are based on Borel, Hilbert transforms, 
positive moments or inverse moments. Sum rule calculations of the decay con- 
stants have been performed since the eighties using Borel transform techniques 
with results within the range f B = 160 - 210 MeV and JbJIb = 1-09 - 1.22, 
|21 EI El HI El HI ■ A more recent sum rule calculations using the new 0{a 2 s ) corre- 
lation function with one heavy and one light current [3] yields higher central val- 
ues f B = 210 ±19 MeV and = 1.16 ±0.04. Lattice QCD determinations 
also give results in a wide range fs = 161 — 218MeV and /b s //b = 1-11 — 1.16 
|lll I12j (for a review and a collection of the results, see A recent HQET 
calculation ^U] yields fs — 206 ± 20 MeV. The large variation of the quoted 
values indicates that there is still room for improvement. 

In general, the sum rule technique assumes duality. In our analysis duality 
means that weighted integrals over experimental measured amplitudes should 
agree with the same integrals evaluated in QCD perturbation theory. 

We use a method based on finite energy sum rules (FESR) which equates 
positive moments of data and QCD theory to evaluate of the decay constants 
of the pseudoscalar bottom mesons {fs and fB,)- The method we propose here 
employs a particular combination of positive moments of FESR in order to op- 
timize the effect of the available experimental information. On the theoretical 
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side it uses the asymptotic (large momentum expansion) of QCD, i.e. an ex- 
pansion in ml/s, where nib is the mass of the bottom quark and s the square 
of the CM energy. Such an expansion makes sense as long as s is far enough 
from the continuum threshold. We will consider the perturbative expansion up 
to second order in the strong coupling constant and up to fourteen powers in 
the expansion of the b-mass over the energy, using the results of the reference 

. On the phenomenological side of the sum rule we will consider the lowest 
lying pseudoscalar _B 9 -meson. In our method we circumvent the problem of the 
unknown continuum data by a judicious use of quark-hadron duality. We use a 
linear combination of finite energy sum rules of positive moments, designed in 
such a way that the contribution of the data in the region above the resonances 
turns out to be practically negligible [T%| . 

The present paper differs from our earlier pilot work |15| in that 0(a s ) 2 
corrections to the perturbative piece and the 0(a s ) corrections to the leading 
non-perturbative term are included. In this way the stability of the prediction of 
fs is improved and the errors are reduced. The stability analysis is also placed 
on more solid footing. 

The plan of this note is the following: in the next section we briefly review 
the theoretical method proposed, in the third section we discuss the theoretical 
and experimental inputs used in the calculation and in the fourth one we write 
up the conclusions. Finally in the appendix we discuss briefly the main issues 
of the polynomials used in the finite energy sum rule. 



2 The method 

The two-point function relevant to our problem is: 

IL(s = q 2 ) = i J dxe^ <n\T(j 5 (x)j 5 (0))\Q>, 

where < | is the physical vacuum and the current j'5 (x) is the divergence of 
the axial- vector current: 

j 5 (x) = (Mq + m q ) : q(x) i 75 Q(x) : 

Mq is the mass of the heavy bottom quark Q(x), and m q stands for the mass 
of the light quarks q(x), up, down or strange. The starting point of our sum 
rules is Cauchy's theorem applied to the two-point correlation function II(s), 
weighted with a polynomial P(s) 

/ P(s) U(s) ds = 0. (1) 

The integration path extends over a circle of radius \s\ = sq, and along both 
sides of the physical cut s G [s p hy S .jSo] where s p h ys . is the physical threshold. 
Neither the polynomial -P(s) nor the power of the integration variable change 
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the analytical properties of IT(s), so that we obtain the following sum rule: 



-/ P(s) ImII(s) = — (f P{s)Tl(s)ds, (2) 

Duality now means that on the left hand side we can use experimental in- 
formation starting from the physical threshold s p h ys . to so, whereas on the right 
hand side, i. e. on a circle of radius sq, we can use the theoretical input given 
by QCD and the operator product expansion. The integration radius Sq has 
to be chosen large enough so that the asymptotic expansion n^ CD (s) of QCD, 
which includes perturbative and non-perturbative terms, constitutes a good ap- 
proximation to the two-point correlator on the circle. 

On the left hand side of equation we can parametrize the absorptive part 
of the two-point function by means of a single resonance B q and the hadronic 
continuum of the bq channel from spoilt. w ith s^^ t s ' ^ s P hys.- Therefore, we 
write the representation of the hadronic experimental data in the form: 

ilmn( S ) = ilmIT cs -( S ) + -ImnX; 9(s - s p J^;), (3) 

7T 7T 7T K J 

where SccSt s t an( is for the physical threshold of the continuum physical region, 
and therefore it is an input in our calculation. It is worth mentioning that it 
must not be confused with the duality parameter (so) that is determined by 
stability requirements in most versions of finite energy QCD sum rules. 

On the right hand side of equation J2J the theoretical input necessary to 
write down the two-point function relevant to the integration over the circle of 
radius so has both perturbative and non-perturbative contributions, 

n QCD ( s ) = n pcrt -( s ) + n nonport -( s ), (4) 

For the perturbative piece, we take the two-point correlation function n pcrt (s) 
which has been calculated in |13| for one massless and one heavy quark as an 
expansion up to second order (three loops) in the strong coupling constant a s 
and as a power series in the pole mass of the heavy quark (M^/s) up to the 
seventh order. We have checked that for the employed range of values of the 
radius sq of the integration contour, the power series converges well and does not 
introduce any appreciable error in the calculation. In the case of one and two 
loops the known complete analytic expressions of the two-point function could 
be used, but we also use the mass expansion in this case because the results of 
the integration can be given the analytically. Numerically it makes no difference 
whether one uses the complete analytic expressions or the expansions. 

The authors of obtain the following compact expansion of the two-point 
function in terms of the pole mass Mb, 

ip^. (s) = n (°)( s ) + ( g^ )n< 1 >( fl ) + ( ^M )V>( S ), (5) 
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where the different terms of the expansion in a s have the form: 
nW( S ),(M i + m,)X £ ( ln 7^y (— V « = M.2)- (6) 



j=-l fe=0 



and the coefficients A^\, of the QCD perturbative expansion are given in ^3] to 
order (M^/s) 7 . For instance, the first term of the expansion is given by: 

n"M-^<iftW.{3-2ta(^) + «^ ta (^ 
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In the asymptotic expansion Q there are also non-perturbative terms due 
to the quark and gluon condensates. In our calculations wc will include these 
terms up to dimension six 



n „o„pc lt . (s) = {Mb + mq f \ Mb (qq) 



(7) 
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The a s correction to the quark condensate [3] turns out to be small but 
non-negligible. 

Before going on, a comment on the asymptotic expansion is in order. It is 
known that the convergence of the perturbative expansion of the two-point func- 
tion, when written in terms of the pole mass, is rather poor. In fact, in many 
calculations involving heavy quarks, the first and second order loop contribu- 
tions are typically of the same order of magnitude making difficult to achieved 
convergence in the results. On the other hand, the expansion in terms of the 
running mass converges much faster over a wide range of the renormalization 
scale. Henceforward, for calculational purposes, we will consider the relations 
among the pole and the running mass in the appropriate order in the coupling 
constant |16M17llT%] in order to rewrite the perturbative piece of order (a s ) 1 (© 
in the form 



Ii^(s) = ml(ii){m b {p)+m q (ii)) 



j=-lk=0 



In - 
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0,1,2) 
(8) 
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and similarly for the non-perturbative piece. The coefficients Aj k depend on 
the mass logarithms ln(TO 2 /^i 2 ) up to the third power. As n(s) is not known to 
all orders in a s , the results of our analysis will depend to some extend on the 
choice of the renormalization point /i. In the sum rule considered here there are 
two obvious choices, \i = mb and /i = sq, the radius of the integration contour. 
The former choice will sum up the mass logs of the form ln(m 2 //i 2 ) and the 
latter choice the ln(— s//i 2 ) terms. We find that the for the choice /j, — mb the 
convergence of the perturbative terms is significantly better, so we will adopt 
this value in the presentation here. 

For the experimental side, we can split up the absorptive part into the es- 
tablished lowest lying pseudoscalar bq resonance B q and an unknown hadronic 
continuum : 

1 lmn(,s) = M% q 4 S(s - M 2 Bq ) + i ImUXs. 9{s - s^) (9) 

where Mb and /b are, respectively, the mass and the decay constants of the 
pseudoscalar meson B q . 

Looking back to equation J2J and taking all the theoretical parameters in- 
cluding the mass of the _B 9 -meson as our inputs in the calculation, we see that 
the decay constant could be computed if we would have accurate information 
on the hadronic continuum contribution, which is, however, not the case. 

To cope with this problem we make use of the freedom of choosing the 
polynomial in equation (J2J. We take for P(s) a polynomial of the form: 

P n (s) = a + ais + a 2 s 2 + a 3 s 3 + . . . + a n s n , (10) 

where the coefficients are fixed by imposing a normalization condition at thresh- 
old 

Pn ( Sph . = M 2 Bq ) = 1, (11) 

and requiring that the polynomial P n {s) should minimize the contribution of 
the continuum in the range sj t s ', sq in a least square sense, i.e., 

s k P n (s) ds = for k = 0,...n-l, (12) 

c phys- 

S CO nt. 

The polynomials obtained in this way are closely related to the Legendre poly- 
nomials. In the appendix the explicit form of the set of polynomials use in this 
work is given. 

This way of introducing the polynomial weight in the sum rule makes negligi- 
ble the integration of i Imfl^*' 6{s — Scont.Oi so that, in the phenomenological 
side of the sum rule only the contribution of the B q resonance remains. On the 
other hand, as we will see in the results, it increases the value of the duality 
parameter sq and, therefore, the asymptotic expansion of QCD can be trusted 
safely in the integration contour. 

To the extend that Im Il^y* - can be approximated by an n degree polyno- 
mial these conditions lead to an exact cancellation of the continuum contribution 



6 



to the left hand side of equation J2J. As a welcome side effect, this choice of 
polynomial will enhance the role of the B q resonance. Notice however that 
increasing the degree of the polynomial P n (s), will require the knowledge of fur- 
ther terms in the mass expansion and in the non-perturbative series. Therefore 
the polynomial degree cannot be chosen arbitrarily high. 

To check the consistency of the method, we have considered a second to fifth 
degree polynomials and the results are compatible within the range of the errors 
introduced by the inputs of the calculation. We also have checked explicitly that 
a smooth continuum contribution had no influence on the result. 

The sum rule method above was previously used in the calculation of the 
charm mass from the cc experimental data. The continuum data in this case 
were known from the BES II collaboration 19s and were shown to have no influ- 
ence on the predicted mass |21j . Employing the same technique, a very accurate 
prediction of the bottom quark mass was also obtained using the experimental 
information of the upsilon system [22J ■ 

After these considerations we proceed with the analytical calculation. The 
integrals that we have to evaluate on the right hand side of the sum rule, equa- 
tion J2J, are 

J {p ,k) = ±J^^) k d s, (13) 

for k = 0, 1, 2, 3 and p = —6, —5, .., n + 1. These integrals can be found e.g. in 
reference [221 ■ After integration, equation (J2J yields the sum rule 

M% q f Bq P(M 2 Bq ) = m2fa)(m 6 (M) + ™ q (p)f (14) 

n 2 6 3 / , n s i 

*EEEE««(^) t&rtfM j( q -j, k ) 

q=0 i=0 j=-l k=0 ^ ' 

+ non-perturbative terms 

where, for brevity, we have not written down the non-perturbative terms explic- 
itly. The contribution of the continuum is neglected. 

Plugging the theoretical and experimental inputs (physical threshold, quarks 
and meson masses, condensates and strong coupling constant) into the sum 
rule, we obtain the decay constant /b. for various values of the degree n of 
the polynomial and various values of sq. Given the correct QCD asymptotic 
correlator and the correct hadronic continuum, the calculation of the decay 
constant should, of course, not depend either on sq or on the degree n of the 
polynomial in the sum rule J2J. Accordingly, for a given n we choose the flattest 
region in the curve /b(so) to extract our prediction for the decay constant. To 
be specific we choose the point of minimal slope. On the other hand, for different 
polynomials, the value of Jb , extracted in this way, could differ from each other 
as the cancellation of the continuum may be incomplete or the QCD expansion 
not accurate enough. We find, however, practically the same results for all our 
polynomials. This additional stability is truly remarkable as the coefficients 
of the polynomials of order n = 2,3,4 and 5 are completely different and the 
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respective predictions are based on completely different superpositions of finite 
energy moment sum rules. This extended consistency leads us to attach great 
confidence in our numbers and associated errors. 

3 Results 

We calculate the decay constants for the B and B s heavy mesons. In the first 
case we take m q = everywhere. In the second case we retain m q — m s in 
the factor (M& + m q ) 2 in front of the correlation function only. Further terms 
in the power series in m 2 /s in (JSJ are completely negligible for the integration 
radii sq we use in the calculation. 

The experimental and theoretical inputs are as follows. The physical thresh- 
old s p h. is the squared mass of the lowest lying resonance in the bq channel. For 
q being the light quark u, we have: 

sph. = Af| = 5.279 2 GeV 2 = 27.87, GeV 2 (15) 

whereas the continuum threshold spoilt * s taken from the next intermediate state 
Btttt in an s-wave I = g, i. e. 

Scont.' = ( m b + 2m,) 2 = 30.90 GeV 2 . 

For q being the strange quark we take: 

s ph . = M 2 Bg = 5.369 2 GeV 2 = 28.83, GeV 2 (16) 

The continuum threshold starts in this case at the value: 

Scant'. = ( m b s + 2m,) 2 = 31.92 GeV 2 . 

In the theoretical side of the sum rule we take the following inputs. The 
strong coupling constant at the scale of the electroweak Z boson mass |2*1] 

a s {M z ) = 0.118 ±0.003 (17) 

is run down to the computation scale using the four loop formulas of reference 
|23| . For the quark and gluon condensates (see for example [Hj) and the mass 
of the strange quark [25] we take: 

< qq > (2 GeV) = (-267 ± 17 MeV) 3 , 
<—GG>= 0.024 ± 0.012 GeV 4 , 

TT 

< qoGq > = ml <qq>, with Wq = 0.8 ± 0.2 GeV, 
m s (2GeV) = 120 ±50 MeV, 

< ss >= (0.8 ±0.3) <qq > . (18) 

As discussed, above we fix the renormalization scale fj, in the theoretical side of 
the sum to be /i = mb(mb), as the perturbative series in a s is well under control, 
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i.e. the first and second order terms are only a few percent of the dominant zero 
order one and the second order one is much smaller than the first order one. 
We use a reasonable variation of fi to estimate the corresponding error in our 
final result. 

Finally, for the bottom quark, a value mb(mb) ~ 4.20 GeV is generally ac- 
cepted. For consistency we take the result mb(mb) = 4.19 ± 0.05 GeV of which 
has been obtained by us [23 with a similar sum rule method. 

In order to calculate the decay constant for the pseudoscalar meson B, we 
proceed in the way described above. We compute fs as a function of Sq with the 
four different sum rules @ corresponding to n = 2, 3, 4, 5. The results, plotted 
in Fig. 1 show remarkable stability properties. We define the optimal value of 
so as the center of the stability region (represented by a cross in Fig.l) where 
the first and/or second derivative of /s(so) vanishes. At these values of sq we 
obtain the following consistent results: 



fs 


= 175 


MeV for 


n = 


2 


fB 


= 177 


MeV for 


n = 


3 


Sb 


= 178 


MeV for 


n = 


4 


Sb 


= 178 


MeV for 


n = 


5 



Notice from Fig. 1 that for the fourth degree polynomial (n=4) there is an 
stability region of about 50 GeV 2 around sq. In this region the decay constant 
changes by less than 1 percent. This stability region is pushed up to 100 GeV 2 
in the case of n — 5. With these considerations we estimate a conservative error 
inherent to the method of ±3 MeV. 

Other sources of errors arising in the calculation of fg are the quark con- 
densates which affect the result by ±3 MeV, and the bottom mass which, in the 
range given above, produces a variation in the decay constant of —13, +10 MeV. 
This last one is the main source of uncertainty in the final result. Finally we 
have changed the renormalization scale in the range /i <E [3,6] GeV. We esti- 
mate an error of ±6 MeV associated to this change in fj, which is roghly related 
to the convergence of the asymptotic expansion 

Other errors due to the QCD side of the sum rule (higher order terms in 
m 2 /s and the error on a s (mz) are negligible in comparison). 

Adding quadratically all the errors, we finally quote the following result for 
the decay constant of the light meson B: 

f B = 178 ± 14 (inp.) ± 3 (meth.) MeV. (19) 

The first error comes from the inputs of the computation and the truncated 
QCD theory whereas the second is due to the method itself. 

Proceeding in the same fashion, but keeping the mass of the strange quark 
in the overall factor and the order m s /mb in the one loop contribution, we find 
the decay constant for the B s meson (/s 3 ), 

f B . = 200 ± 14 (inp.) ± 3 (meth.) MeV. 



9 




Figure 1: Decay constant fg as a function of the integration radius so for 
mb(mb) = 4.19 GeV. The lines represent the sum rule for several choices of 
the polynomial P„(s). 



(compare Fig. 2) 

In the analysis of theoretical errors the only new ingredient is the uncertainty 
coming from the strange quark mass which turns out to be negligible. 

The ratio of the decay constants fg B and fg (which would be 1 in the chiral 
limit) is of special interest. We find: 

^ = 1.12 ± 0.01 (inp.) ± 0.03 (mcth.). (20) 

JB 

In the calculation, the uncertainties due to the method and to the theoretical 
inputs (mainly to the bottom quark mass) are correlated, so that the final error 
is very small. 

Finally we should mention that in the past the stability region was often 
determined in a rather ad hoc fashion by considering the intercept of sum rule 
predictions for moments differing by one power. We find that, suitably modified, 
this prescription also works in our sum rule method with similar results but 
larger errors. 
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Figure 2: Decay constant /b s as a function of the integration radius s for 
mb(m,b) = 4.19 GeV. The lines represent the sum rule {2J for several choices of 
the polynomial P n (s). 



4 Conclusions 

In this note we have computed the decay constant of -B 9 -mesons for q either 
the strange or the u or d massless quarks. We have employed a judicious com- 
bination of moments in QCD finite energy sum rules in order to minimize the 
shortcomings of the available experimental data. On the theoretical side of 
the pseudoscalar two-point function, we have used in the perturbative piece an 
expansion up to three loops in the strong coupling constant and up to order 
(ml I s) in the mass expansion and in the non-perturbative piece we consid- 
ered condensates up to dimension six including the a s correction in the leading 
term. Instead of the commonly adopted pole mass of the bottom quark, we use 
the running mass to get good convergence of the perturbative series. It turns 
out that for the renormalization point fi = mbirnb) the first and second order 
contribution in the strong coupling are term by term much smaller. This good 
convergence is due to the summing up of the mass logarithms. 

In the sum rule, the contour integration of the asymptotic part is performed 
analytically. This particular fact differs from other computations based on sum 
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rules where the asymptotic QCD is integrated along a cut of the two-point func- 
tion starting at the pole mass squared. The latter way to proceed is problematic 
when loop corrections are included and the complete analytical QCD expres- 
sion along the cut is not known. In this approach QCD has to be extrapolated 
from low energy to high energy [3]. We also differ from many other sum rule 
calculation of fs in that we do not require two largely unrelated sum rules to 
determine a stability point via an intercept of the curves /s(so)- 

Our results are very sensitive to the value of the running mass, giving most 
of the theoretical uncertainty. On the other hand they turn out to be very stable 
against the variations of the other parameters, in particular the renormalization 
scale and the integration radius sq. 

Comparing with other results in the literature, our results agree within the 
error bars with the ones obtained using sum rule methods and with lattice 
computations. When compared with the most recent numbers pfl 1111 IT2l IIP) , 
however, our results are lower. 

The interest in our method lies in the fact that it approaches the problem 
from a different angle and, in our opinion, is less proned to systematic errors. 



Appendix 

For convenience of the reader in this appendix we list the first few polynomials 
emerging from relations i|llll2|) . From the second condition, namely l|12l) . it 
is easy to realize that the set of polynomials P n (s) are n-degree orthogonal 
polynomials in the interval of the variable s £ [scoiit 8 j s o] ■ Then, leaving aside 
the normalization condition (|ll|l that we take, for convenience, in order to stress 
the contribution of the lowest lying resonance in the sum rule, they are related 
to the so-called Legendre polynomials (P n (x)) in the interval of the variable 
x G [—1,1]. After including this normalization condition, we can write more 
precisely: 

Pn(s) = V ? iX{S l\ (2D 

Where the variable x(s) is: 



x(s) = 



phys. 
s — s cont. 



which ranges in the required interval when s £ [Scont. > s o] ■ 

The explicit form of these polynomials is well known and can be found, for 
instance, in [2EJ. Nevertheless, for sake of completeness, we quote here the ones 
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Figure 3: Polynomials P n (s) obtained from conditions I|llll2fl taken at the 
stability values of Sq in the calculation of fs- 

we have used in the calculation. 

V 2 (x(s)) = \{x*-l), 

-Ps(x(s)) = i(5z 3 -3x), 

Va{x{s)) = i(35a; 4 - 30x 2 + 3), 
8 

V 5 (x(s)) = l(63x 5 - 70x 3 + 15a;). 
8 

Finally, in Fig. 3, and in order to appreciate how the dumping of the ex- 
perimental physical continuum in the sum rule is expected, we have plotted the 
form of the polynomials P n {s) for n=2,3,4,5 at the stability values of sq used in 
the calculation of Jb- 
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